Predicting spinor condensate dynamics from simple principles 
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We study the spin dynamics of quasi-one-dimensional F — 1 condensates both at zero and finite 
temperatures for arbitrary initial spin configurations. The rich dynamical evolution exhibited by 
these non-linear systems is explained by surprisingly simple principles: minimization of energy at 
zero temperature, and maximization of entropy at high temperature. Our analytical results for the 
homogeneous case are corroborated by numerical simulations for confined condensates in a wide 
variety of initial conditions. These predictions compare qualitatively well with recent experimental 
observations and can, therefore, serve as a guidance for on-going experiments. 

PACS numbers: 03.75.Mn, 03.75.Kk 



Spinor condensates realized by optically trapped ultra- 
cold atoms in a hyperfine Zeeman manifold allow to ad- 
dress a broad scope of problems, which are related to 
magnetic ordering. These include quantum phase tran- 
sitions, exotic topological defects and spin domains ei- 
ther within a mean-field regime or within the strongly- 
correlated one [ll. Already in the mean-field regime the 
dynamics of spinor condensates shows some intriguing 
features, with an apparent randomness regarding the evo- 
lution from a given initial state towards a steady state, 
accompanied by formation of spin domains. Even for the 
simplest case of an = 1 spinor condensate, the inter- 
play between spin-spin interactions, non-linear terms and 
temperature effects rends the analysis of the dynamics 
rather complex. Previous studies of the spinor dynam- 
ics show that, at early stages of the evolution, there is 
a coherent population transfer between the different hy- 
perfine coupled sublevels 0, Q ■ Inclusion of temperature 
(T) not only smears out the population transfer Q , but 
also leads to a different distribution of population among 
the different hyperfine levels. 

In this contribution we analyze the dynamics of spinor 
_F = 1 condensates both at zero and finite T from a new 
perspective. As we shall show, the complex dynamics 
displayed by these systems can be understood in terms 
of oscillations in phase space around a steady state. The 
configuration of this state can be approximately deter- 
mined by analyzing the trajectories of constant energy 
in the homogeneous case. To a good approximation, the 
populations that characterize this state are rather close 
to those that minimize the energy associated to spin-spin 
interactions for a given magnetization. In contrast, at fi- 
nite T, the system is able to exchange energy with the 
thermal clouds. At large enough temperature, the popu- 
lations of the steady state can be simply determined by 
maximizing the entropy of the homogeneous condensate. 
This leads to a different long-time configuration of the 



populations than the one attained starting from the same 
initial conditions at T = 0. Our claims are supported by 
analytical results for an homogeneous condensate and by 
numerical investigations for trapped condensates. They 
provide a good qualitative agreement with recent experi- 
mental data on the dynamics of spinor condensates [j, 0] . 
Furthermore, our results can be straightforwardly used to 
analyze and predict experimental outcomes. 

In the mean-field approach, the spinor condensate is 
described by a vector order parameter \E' whose compo- 
nents ipm correspond to the order parameter of each mag- 
netic sublevel \F — I,m) with m ~ 1,0,-1. In absence 
of an external magnetic field and at T = 0, the properties 
of the condensate are determined by the spin-dependent 
energy functional 0, @, 01 : 



where repeated indices are assumed to be summed. The 
trapping potential is taken to be harmonic, axially- 
symmetric and spin independent. The coefficients cq and 
C2 describe the spin-independent and the spin-dependent 
contributions of the binary elastic collisions of two iden- 
tical spin-1 bosons. They are expressed in terms of the 
s-wave scattering lengths a/ of the combined symmet- 
ric channels of total spin / = 0,2 as Cq = 47r?i^(oo + 
2a2)/(3M) and C2 = ATTh'^{a2 - ao)/(3M), where M is 
the atomic mass. F are the spin-1 matrices 0,0] ■ In this 
Letter, we focus on ^''Rb spinor condensates 0,0, [1,0], 
which have ferromagnetic character (c2 < 0). 

According to ihd^m/dt — S£/S^p*^, a set of coupled 
dynamical equations for the spin components are ob- 
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tained: 



ihdil)±i/dt = \Hs + C2(n±i + no - n=pi)] i>±i 



(2a) 



i-i)] 



(2b) 



+ C2 ^AoV-i^l > 

ihdipo/dt = [Tis + C2(ni 4 

+ 2c2?AlV'oV^- 

with = ~h^/{2M) + Vext + Con being the spin- 
independent part of the Hamiltonian. The density of 
the m-th component is given by |'0m(r)P, while n(r) = 
X^m IV-'m(r)P is the total density normahzed to the total 
number of atoms N . The population of each hyperfine 
state is N„i ~ J dr\ip„i{r)\'^ . Defining the relative pop- 
ulations A„i — N„i/N, the magnetization of the system 
is Al = Ai — A_i. The total number of atoms and the 
magnetization are both conserved quantities [2]. Then, 
the conditions Ai -I- Ao -I- A_i = 1 and A4 — Xi — A_i, 
together with the positiveness of Am, restrict the allowed 
configurations of the system in the {Xo,A4) phase space 
to the region below the dotted line in Fig. [TJ 

To understand the physics governing the spin dynam- 
ics, it is instructive to analyze first the homogeneous sys- 
tem at T = 0. Translational invariance leads to the fol- 
lowing ansatz: ip„i = V^V Xm e'^™ . Notice that, since 
the system is homogeneous, the value of n is not rele- 
vant to determine the relative populations that minimize 
the energy. Thus, the only relevant contributions to the 
energy are those associated to C2'. 

E,2iXo, M, 9) ^ C2n^ [Xo{l ~ Xo) + My 2 



+ Ao v/(l - Ao)2 - cos e], (3) 

where 9 = 29o — 9i — 9-i. For a given value of A4, 
the ground state configuration is obtained by minimiz- 
ing Eq. ([3]) with respect to Aq and 9. For ferromagnetic 
systems, this yields 9 — and 

k^-M'), (4) 



which defines the solid line drawn in Fig.[TJ For instance, 
in the symmetric case {M = 0) the minimization leads 
to a ground state with Aq = 0.5 and 9 = 0. 

We study now the dynamical evolution of the homoge- 
neous system by solving numerically the set of Eqs. 
after removal of the terms —h"^/ (2M) + Vext- We con- 
sider three independent phases 9m in the time evolution, 
although the energy depends only on the relative phase 
9 [Eq. ([3])]. As expected, we find that starting from an 
equilibrium configuration {e.g. Xq — 0.5 and 9 = for 
= 0) there is no evolution at all. However, if the 
initial configuration does not correspond to the equilib- 
rium state, then the system evolves exchanging popula- 
tions between spin components, but conserving energy 
and magnetization. 

Therefore, in the homogeneous system the dynamical 
evolution of the spin populations is characterized by tra- 
jectories of constant energy and magnetization in the 
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FIG. 1: (color online) Equilibrium configurations in the 
Ao — plane for the homogeneous condensate at T = 
[solid line, Eq. (|4])] and at large T [dot-dashed line, Eq. ((5])]. 
The dashed line in between indicates Aq"'", the geometrical 
center of the closed orbit with the same magnetization, while 
the dotted line bounds the region of physically allowed values 
for (Ao, M). Full symbols stand for the time averages of Xo{t) 
calculated from the numerical time evolution for different ini- 
tial configurations at zero and finite T for the inhomogeneous 
system; open squares indicate the initial conditions for the 
simulations ai M = 0.2, 0.4 respectively. Due to the sym- 
metry of the equations, it suffices to present the results for 
positive magnetizations. 



Xq — 9 plane, which are defined by the starting values 
of A^,Ao and 9. In most cases, i.e. if the initial condi- 
tions are not very far from the equilibrium configuration, 
defined by Ap' and 9 = 0, the trajectories in the Aq — 6* 
plane are closed orbits (see also [11]). Then, the evolu- 
tion of the populations -in particular of Ao- shows up as 
oscillations of a well defined frequency, around the ge- 
ometrical center of the orbit. An analysis of the orbits 
defined by Eq. ([3]), shows that these centers lie always be- 
tween Ao' and AJ,""' = (1 - |Al|)/(2 - |7W|), represented 
by the dashed curve in Fig. [TJ With this information in 
mind we move now to the confined systems. 

To this aim, we solve the equations of motion for 
20000 atoms of *^Rb, in a highly elongated trap llJ± = 
2n X 891 Hz and cj^ = 2Tr x 21 Hz, i.e. wj_ '> uj^, 
so that Eqs. ^ become quasi-one-dimensional. The 
corresponding scattering lengths are ao = 101. Sas and 
02 = 100. 4aB, and the coupling constants co and C2 
are accordingly rescaled by a factor l/{2TTa']J, being 
a_L = y/h/mu>± the transverse osciUator length [loj . 
The validity of this approximation for the case of a spin-1 
atomic condensate has been proved for a very wide range 
of conditions in Ref. ll|. We take as initial wave func- 
tions 'ipmiz, t = 0) the wave function of the ground state 



of the scalar condensate, obtained from Tig with coupling 
constant co/(27ra5_), and normalized to their correspond- 
ing initial populations Nm- This corresponds to a Single 
Mode Approximation (SMA) ansatz for the initial condi- 
tions. SMA provides a good description of spin systems 
in certain conditions [l3 . [isj . However, in solving the 
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time evolution equations we keep the full dependence of 
on position and therefore we go beyond the SMA. 
Actually, this is a necessary condition to study the for- 
mation of spin domains. 

The numerical solutions of the dynamical equations 
show oscillations of the integrated spin populations 
around a steady state, which coincides with that in the 
homogeneous case [Eq. ([4])]. Indeed, if we start with 
a configuration lying on the curve defined by Aq' with 
= 0, the system does not evolve at all. However, 
any initial configuration with (Ao,A^) on the curve, but 
6* 7^ 0, yields dynamical evolution of the spin popula- 
tions around the corresponding steady configuration of 
the homogeneous case. 

Our results are summarized in Fig. [T] The dots cor- 
respond to the time average value of \o{t) obtained by 
solving the equations of motion for different initial out- 
of-equilibrium values of A^, Aq, and 9 = Q. All cases 
studied would correspond in the homogeneous system to 
closed orbits in the Xq — 9 plane. The magnetization is 
constant throughout the evolution, while Xoit) presents 
slightly damped oscillations around the same Aq as in 
the homogeneous case, which is restricted to be in the 
narrow region limited by the solid and the dashed curves 
described above. Thus, the curve Ag'' provides a reli- 
able estimation of the configuration of the steady state. 
Notice that in all reported cases, the dots lie very close 
to the curve X'^ as we have taken, in the selection of the 
starting conditions for ^ 0, configurations and phases 
not far from the equilibrium conditions. 

As an example of the dynamical evolution of the sys- 
tem at T = 0, we display in Fig. [Ifa) the spin dynamics 
for an initial symmetric spin configuration correspond- 
ing to Ai — A_i — 5% and 6* = 0. The oscillations 
between m = and m = ±1 populations are not regular 
and present a dynamical instability around tdom, when 
the condensate starts the multidomain formation process 
[3, [3, 3, l3| • An estimate of the time scale for the ap- 
pearance of the instability is tdom ~ 27rS/(c2no) (isl. 
where no is the density at the center of the trap. In our 
simulation, no ~ 2 x 10^^ cm~'^, leading to tdom ^ 140 
ms, which is in agreement with the dynamics of the sys- 
tem. Around this time, the condensate starts the for- 
mation of small dynamical spatial spin structures, also 
observed experimentally [1, [J, while the relative popula- 
tions oscillate around the equilibrium state (Aq'' = 50%, 
A^j = 25%). Due to this inhomogeneity of the phases 
and density profiles , it is not possible to associate the 
dynamics to a well-defined orbit in the Xq — 9 plane. 

It is interesting to note the formation of spin domains 
despite the total density profile remaining almost un- 
changed [3, . Hence, spin domains are hnked to the 
kinetic and spin-spin energy terms, since all the other 
contributions to the energy functional depend only on 
the total density. The evolution of the kinetic energy 
and Ec2 is shown in Fig.[2Kc) as a function of time for the 
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FIG. 2: (color online) Spin dynamics for the initial configura- 
tion (Ai, Ao, A-i) = (5%, 90%, 5%) and 6* = 0. Relative popu- 
lations as a function of time for (a) T = and (b) T = 0.2rA • 
(c) Evolution of the kinetic (black) and Ec2 (red or dark grey) 
energies, and their sum (green or gray) at T = 0. The vertical 
lines are at tdom- 



first 500 ms. Clearly the appearence of domains is linked 
to a rapid variation of kinetic and E(.2 energies. More- 
over, since the system has no dissipation at T = 0, the 
total energy has to be conserved, and therefore, the en- 
ergy difference between the starting and the equilibrium 
configurations provides an upper bound to the energy 
associated with the formation of spin domains. Actually 
this energy is always rather small, as can be seen in the 
figure. In the case of Fig. [S] the excess energy amounts 
to Q.2^hhJz which represents a 0.3% of the total energy. 

We incorporate now finite-temperature effects by 
means of a Bogoliubov-de Gennes description of the ther- 
mal clouds. Further we assume that each spin component 
has its own cloud [l7| , and solve the corresponding time- 
dependent equations under these conditions . With our 
choice of parameters, the condensation in an elongated 
trap occurs at ksT^^ « O.Q^hujhoN^^^ ^ where Who = 



(cj^cj^)i/3 [18]. In Fig. EKb) we show for T = 0.2TA 
the time evolution of the relative populations, starting 
from the same initial parameters as in the T = case. 
It is worth noticing that even if T > hw^, the quasi- 
onc-dimcnsional description of the thermal effects is still 
vahd |.19.] . 

A simple comparison between Fig. 2(a) and 2(b) shows 
that at finite T the oscillations in the populations are 
clearly damped and that the steady state configuration 
that is reached is different from the one at T = 0. In 
particular, in Fig.l^jjb) we observe that starting from the 
same initial configuration (5%, 90%, 5%), the long-time 
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population distribution changes from (25%, 50%, 25%) at 
T = to a steady configuration with all m-states al- 
most equally populated (equipartition). Moreover, even 
though the local magnetization is no longer conserved as 
it was at T = y, the total magnetization is still a 
conserved quantity along all the time evolution. In gen- 
eral, temperature effects reflect in a faster damping and 
in an earlier domain formation in comparison to T = 0. 
A way to understand this result is to assume that the 
steady configuration is attained by maximizing the en- 
tropy, subjected to the conservation of total population 
and magnetization: 

5 = -^A™ ln(A„0-A(^mA„-A^)-/x(^A™-l), 

■m m m 

where we have introduced the Lagrange multipliers A and 
H associated with the magnetization and the normaliza- 
tion, respectively. Imposing that dS/dXm = 0, together 
with the normalization and magnetization constraints, 
leads to the following condition for the steady state con- 
figurations: 

= l-2Xa- 3X1, (5) 

which corresponds to the dot-dashed curve plotted in 
Fig. [T] Notice that the curve thus obtained is indepen- 
dent of temperature. However, the hypothesis underly- 
ing its derivation is appropriate for a high temperature 
regime. Therefore, one expects that this curve together 
with the one for the T = case [Eq. ([4])] limit the zone of 
the configurations around which the system will oscillate 
depending on temperature. In fact, the results of our 
numerical simulations for different temperatures and ini- 
tial configurations, shown in Fig. [l] are spread between 
the two curves. We present results for T/T^^ = 0.05 
(rhombi), 0.1 (squares), 0.2 (triangles) and 0.4 (inverted 
triangles). For a given initial configuration, when the 
temperature increases, the steady state gets closer to 
the curve of maximum entropy. This is illustrated for 
A4 = 0.2 and 0.4. In both cases, the initial config- 
uration is indicated by open squares. Notice that for 
A4 = 0.2, the initial configuration is on the equilibrium 
curve (T = 0) and its evolution is a direct consequence of 
the temperature. Also one observes a larger dependence 
of the final steady state configuration on the initial con- 
ditions than at T = 0. In the figure, this is illustrated for 
A4 — 0, and 0.1 by the presence of more than one sym- 
bol belonging to the same temperature that correspond 
to different initial conditions. 

Summarizing, we have shown that at zero temperature 
and in absence of external magnetic fields the spin pop- 
ulations of = 1 condensates oscillate around a steady 
state which is well determined by analyzing the simpler 
homogeneous case. This steady state is close to the one 
described by minimizing the Ec2 energy at constant mag- 
netization for the homogeneous case [Eq. Q]. In confined 



systems, spin domain formation is associated to the ex- 
cess of spin-spin coupling energy, which converts into ki- 
netic energy back and forward dynamically. 

At finite T, the condensate exchanges energy with the 
thermal clouds, and the system evolves towards a differ- 
ent steady state configuration which depends also on the 
temperature. We have shown, however, that the configu- 
rations defined on the homogenous case by minimization 
of the energy and maximization of the entropy set an 
upper and lower bound on the accessible steady state 
configurations for confined spinor condensates at finite 
T. These predictions can serve as a guide for the inter- 
pretation of on-going experiments on spinor condensates. 
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